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Two non-interacting quantum systems which couple to a common environment with many degrees 
of freedom initially in thermal equilibrium can become entangled due to the indirect interaction 
mediated through this heat bath. I examine here the dynamics of reservoir induced entanglement 
for a heat bath consisting of a thermal electro-magnetic radiation field, such as black body radiation 
or the cosmic microwave background, and show how the effect can be understood as result of an 
effective induced interaction. 
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I. INTRODUCTION 

About two decades of research in quantum information have led to the picture of quantum entanglement as a precious 
resource that plays a key role in processing information more securely and more efficiently than classically possible LJ 
Entanglement is necessary for quantum teleportation and the exponential acceleration of quantum algorithms 
, and it allows for secure quantum key distribution |4| . Recent experimental demonstrations of quantum teleportation [_ 
IH3 an d small scale quantum computation 0] confirm this picture. Creating and manipulating entanglement 
("O \ in a controlled way remains a challenge, as environmentally induced decoherence tends to rapidly destroy entanglement. 
00 ■ On the other hand, it has been shown that the coupling of two quantum systems to a common heat bath can also 
create substantial (mixed state) entanglement. This effect was originally demonstrated in the framework of an exactly 
solvable model |l2j. and has been confirmed by perturbative calculations in the Markovian regime 0,^3- Cirone et 
al. have unveiled a connection to the Casimir-Polder interaction |15| , and Oh and Kim have shown by a renormalization 
group analysis how the effect can be understood through an effective induced interaction between the quantum systems 
via the common heat bath [l6l | . 
I Heat baths are ubiquitous, and one might therefore wonder, if "reservoir induced entanglement" (RIE) |17| . i.e. cre- 
ation of entanglement through coupling to a common heatbath is common place and over what distances and with 
what time dependence it might arise. For example, all quantum systems containing charged particles couple to the 
ambient thermal electro-magnetic radiation, i.e. the black body radiation (BBR) in a laboratory or cryostat at the 
corresponding temperatures, and even in free space there is still the cosmic micro-wave background (CMB), a basi- 
cally perfect black body radiation at an absolute temperature T = 2.728 ±0.004 K that fills the entire known universe 
[21). Do these heat baths induce entanglement between remote quantum systems? 
• t-H ', In this paper I show that BBR effectively constitutes two different heat baths which couple differently to the couple 
Alice-Bob, and whose effects largely cancel when it comes to entanglement creation. Entanglement is therefore 
created only very slowly, far behind the light-cone, and the entanglement created oscillates as a function of time. The 
entanglement can be close to perfect, but the first maximum entanglement arrives only after a time t\ which scales 
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like the 3rd power of the distance R between the two quantum systems, 
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where ao ~ 1/137 and cq are the fine structure constant and the speed of light in vacuum, respectively, and d denotes 
the dipole moment of the quantum system divided by the electron charge. The slow creation of entanglement limits 
the distance over which it can be created before competing decoherence processes set in. 

I will elucidate the role of the high-frequency, far off-resonant modes of the heat bath, and discuss the temperature 
dependence of the phenomena. 



II. THEORETICAL FRAMEWORK AND PHYSICAL SYSTEMS 



Let us consider two identical quantum systems A and B which couple to the thermal electromagnetic radiation 
field in open space, as exemplified by the CMB, but which do not interact directly. We assume that A and B can be 
approximated as two state systems with states |0) and It turns out that the standard quantum optics approach of 
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rotating wave approximation and Markovian Master equations based on Fermi Golden Rule rates (i.e. second order 
perturbation theory in the atom-field coupling constants) is not fully adequate for describing RIE. Firstly, the explicit 
dependence of the time scale on which RIE is produced on the cut-off frequency of the heat bath found in ^3] hints 
to the importance of non-resonant modes, which invalidates the rotating wave approximation. One might argue that 
in [12| non-resonant modes came into play because degenerate energy levels were considered, but even for finite level 
spacing A the high-fre quen cy modes should be relevant for times t -C A. Note that the rotating wave approximation 
was also avoided in (lUllSj. 

Secondly, there is evidence that for the specific heat bath and quantum systems to be discussed here RIE is an effect 
that arises only at fourth order in the coupling constants • Both issues combined call for a fourth order calculation 
without rotating wave approximation, which makes the theory very heavy. On the other hand, the problem can be 
solved exactly and with less effort in the case of degenerate energy levels [l3.l23| . This approach will be followed here. 

From an experimental point of view, two-state systems with exactly degenerate energy levels are hard to find. 
However, it should be kept in mind that a finite level spacing just introduces an upper limit to the time for which 
the present theoretical analysis is applicable: for any experiment terminated within a time i « 1/A the system 
hamiltonian of the two-state systems A and B can be neglected [2]|. Whether or not entanglement can still be 
produced beyond this time is an interesting experimental question. 

Double quantum dots (DQDs) seem to be a promising candidate, and the following analysis will be geared specifically 
towards these systems. Recently, the coherent manipulation of two states |0) and |1) located in the two wells of a 
single such device, as well as state preparation and state measurement were demonstrated; coherence times of the 
order I-10ns were achieved The energy barrier between the two wells as well as the energies of |0) and |1) in 

a DQD can be tuned with the help of gate voltages and the barrier can be made very high after the preparation of a 
superposition, such that |0) and |1) become to good approximation degenerate eigenstates of the DQD hamiltonians. 
The question of whether or not the two state approximation still holds when we have to consider high-frequency 
modes might be addressed eventually experimentally by trying to shield each DQD with its own superconducting 
cavity. Modes with frequency higher than the superconducting gap will be absorbed and one might thus envisage to 
control the cut-off frequency of the heat bath. Otherwise a cut-off of the order Tilo ~ leV arises naturally due to the 
band-gap of the semi-conductor material in which the DQDs are embedded; e.m. waves with higher frequencies get 
absorbed in the semi-conductor. A disadvantage of the DQDs are the competing intrinsic decoherence mechanisms 
such as phonon scattering [2(| and fluctuating electric fields j2?J other than those from the BBR, which will be 
neglected in the following analysis. 

As for the heat bath, I will consider specifically BBR with periodic boundary conditions. One might wonder if the 
mode structure (and thus the boundary conditions) make a difference, as each mode couples to two spatially separated 
quantum systems. For the geometrical situation considered below, and within dipole approximation of the coupling, 
it turns out that a box-shaped cavity with perfectly conducting walls leads, in the limit of infinite volume and fixed 
distance R, to the same interaction hamiltonian, and one can thus read in the following CMB with periodic boundary 
conditions or BBR in a box-shaped cavity interchangeably. 



III. DOUBLE QUANTUM DOTS INTERACTING WITH BBR 

Let us arrange the two DQDs such that the axes of the dots are aligned with the vector joining them, designated 
as z axis in the following. The position operators of an electron in dot 1 and 2 have the matrix elements (0|zi ; 2|0) = 
—d/2 = — (1 1 ^1,2 1 1} ■ All other matrix elements of ,21,2 vanish to very good approximation due to the exponentially 
small overlap of the states |0) and |I), and so do the matrix elements of the other electron coordinates Xi$ and j/i^ 
for an assumed even parity of the ground state wave functions. 

For describing the BBR I use a slightly unconventional representation of the electro-magnetic field, which is very 
useful for the exact treatment for vanishing level spacing (or in general, if one does not use the rotating wave 
approximation). In fact, the BBR can be considered as two independent heat baths (see appendix lVIII|) . one containing 
cos-waves, the other sin-waves, with an electric field operator 



E (M) = ^2 ^fce k a(QkaiCoskr + Q ka2 sinkr) , (2) 

and canonical coordinate operators Q\^ al/ of the harmonic oscillators corresponding to the different electro-magnetic 
field modes. The k are quantized wave vectors {ki = 27rrij/i, n, e Z for i = x, y, z, for quantization in a volume 
L 3 with periodic boundary conditions), a = 1, 2 counts polarizations, v distinguishes the cos-waves {v = 1) from the 
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sin-waves (y = 2), and the ' on the sum means summation restricted to k x > 0. The free field hamiltonian reads 

I 1 

-ffbath = X! X! ( P kav + u lQkau) > ( 3 ) 

k a=l,2w=l,2 

with ojk = co|k| = coA;. 

Placing the two DQDs at positions R/2 and — R/2 (R = (0, 0, R)) we obtain in dipole approximation the coupling 
hamiltonian 

Hmt = 2^ (( <Jzl + a zz) ffkalQkctl + (o"zl ~ <Jz%) 9\ialQ\iOLl) , (4) 

k a=l,2 

written in terms of Pauli matrices a z i, a z2 in the basis (|0), |1)) of dots 1 and 2. The coupling coefficients g-^av are 
given by 



5kai = -pr\ /tt— w fe e kQ • cos — , (5) 




and the same equation holds for (?ka2 up to the change cos — ► sin. 

In the regime of energy degenerate states |0) and |1) discussed above, the time evolution of each dot due to its 
own system hamiltonian can be neglected, and the total hamiltonian is thus simply H — H lnt + -ffbath 12, 23]. The 
interaction hamiltonian represents a generalization of the situation considered in |23(, in the sense that each mode 
of the BBR couples through a different "system coupling agent" Skew = gkau(&zi — (— l)"^)- Note, however, that 
all coupling agents commute with each other. This allows to generalize the time evolution derived in 0,H^| of the 
reduced density matrix describing the two DQDs alone to 



i\p(t,t )\s') = expj-^ f (A S)kai / - A S ',kai/) 2 f(k) (6) 
-i - A^kaJ ¥>(*)] }(s\p(0, t )\s'), 



where \s) is one of the four states 1 00) , 1 01) , |10), or |11), and A Sj kcw are the corresponding eigenvalues of Skew- It was 
assumed that the bath is initially in thermal equilibrium, and that the total initial density matrix factorizes into a 
system part and a bath part. The dependence on the time of travel of a light signal between the two DQDs, to = R/cq 
with R = |R|, will appear below. 

For the arrangement of the DQDs described above, only one polarization direction (a = 1) contributes, with 

£ki = £0 in polar coordinates for k. With the abbreviation Gk = — ^ j^-uJk sin 6 the relevant eigenvalues are 

Aoo.kii = 2 cos(fci? cosd/2)Gk = — An )k ii and Aoi, k i2 = 2sin(fci?cos#/2)Gfc = — Aio,ki2- All other eigenvalues vanish. 
The functions f(k) and ip(k) depend on k = |k| as (fi = l/^^T)) 

(3hujk\ 1 — cosajfct 1 / sinwfct 



m = coth iv ) -mjr> ^rvj- (7) 

Transforming the sums over modes k into integrals for the limit of large L, we find 

( Sl \p(t,t Q )\s 2 ) = exp ^- A(/i(i,io)C ai , a2 + h(t,t )S SuS2 -i((pi(t,to)C auaa + tp 2 (t,t )S sllS2 ))j{ Sl \p(0,t )\s 2 }^) 
with 
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S -\^An^\l G -\^na^\^ S -\^(\ n 1 > C ~ ~ S ■> ( 9 ) 



in the basis |00), 1 01) , |10), |11), and A = a \ 2 , where ao = e 2 /4ire Tic ~ 1/137 is the fine-structure constant. 

irc o T 

In the following, both t and to will be expressed in units of the thermal time r = (3h for all finite temperatures. The 
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functions /i, 2(^,^0) and f 1.2(^,^0) are then given by 

MM = £ _ "*,,coth( s /2)(l-co aW )(I + (-ir(^-2g|i)) (10) 

The integrals extend in principle from zero to infinity, but a UV cut-off y max = w max T is needed to regularize them. 

Due to C = —S, the final density matrix depends only on the phase difference tp-(t) — (fi(t) — if2(t), which can be 
written in closed form as 

<f-(t,t ) = -3 (-2 sin {y max t Q ) + Si {y max {t - h)) + 2Si {Vmaxh) - Si(y max (t + t Q ))) 
2 

H -5 sin(y ma xt)sm(y max to) ■ (12) 

ymax^Q 

Note that (i, t ) remains finite for y max — > 00. For large t and to, ip-(t,to) increases proportional to t and decays 
as l/£g. The functions f v {t,to) on the other hand scale like y^ax with the cut-off, which gives a physical significance 
to the high-frequency modes. The influence of the form of the cut-off function will be examined below. For the 
moment we assume the simplest form, a sharp cut-off at u> = io max . 



IV. RESERVOIR INDUCED ENTANGLEMENT 



I have evaluated the entanglement of formation E{p(t,to)) for an initial state in the form of a pure product state, 
(|0) + |1))® (|0) + |l))/2, i.e. E(p(t,t )) — using Wootter's formula [2^ and by numerically integrating the remaining 
expressions for f\^{t,to). Alternatively, one may approximate coth ~ 1 for y ma x S> 1, which allows for an analytical 
solution also for /„. The result for T — 2.73K and d = lOnm is shown in Fig^ Clearly, there is no entanglement for 
space-like separated points, t < to. However, almost perfect mixed state entanglement arises for t > to, starting at 
around t ~ 100r. The entanglement oscillates rapidly as function of t, and vanishes for t/r < c(to/r) 3 , where c is a 
constant of order 10 12 in agreement with eq.Q. This behavior can be understood by considering the long time limit 
i > to > t, which leads to /1.2 — y^ax/® f° r Umax 3* 1- In this case the final density matrix depends only on two 
variables, v = duj max /co and 0_ — Aip- — </>_(t,io), with 

( S \p(v,^)\s') ~ exp(-^-v 2 (S s , s ,+C s , s ,)+i^C S:S/ ) (s\p(v,<P-(0,t o ))\s') . (13) 

Fig0 shows the corresponding entanglement. The function is 7r-periodic in 0_ with a first maximum at 0_ = tt/2, 
which leads to eq.Q for the first maximum entanglement. The scaling of t\ with i? 3 makes the entanglement 
production extremely slow: during the lifetime of the universe, it would have reached a distance of only about 8.4km 
for d = Ipm. For the same dipole moment, a distance of about 52pm should be reached during a coherence time of 
the order of 100ns. The maximum amount of entanglement as well as the time to first finite entanglement creation 
is controlled by v, which has to be smaller than about 50. One should therefore try to have a large d for a large 
maximum distance and a small to rnax to get large maximum entanglement, whereas the temperature becomes irrelevant 
in this regime (huj max 3> k^T). Note that r cancels in the expression for for large t, to, such that also t\ becomes 
independent of temperature. 



V. THE EFFECTIVE INTERACTION 



The physical origin of the entanglement generation can be understood as arising from an effective interaction H c ff 
mediated through the coupling to the common bath |16(. This idea can be made quantitative by observing that a 
Hamiltonian of the type 



k k 



(14) 
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FIG. 1: Entanglement of formation E for two initially not entangled DQDs with d — lOnm coupled to the CMB at T — 2.73K 
as a function of log 10 (fo/r) and log 10 (t/r), r = j3Ti. Black means perfect entanglement, E = 1, white no entanglement E — 0. 
Entanglement is created only for t/r ^ 10 12 (to/T) 3 - 
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50 



FIG. 2: Entanglement of formation for t > to > r as function of v = uimaxd/co and the phase <p-. The latter increases linearly 
with time t leading to entanglement oscillations, and decays as t^ 3 = (R/co)~ 3 , which makes the entanglement creation very 
slow for large distances R. 
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where S k is an arbitrary operator ac ting on the Hilbert space of the system (i.e. here the qubits of Alice and Bob) 
leads to the time evolution operator |29| 

e -i/HHt = D{ _ £ Et s t k)e ^MH b , tM!f )t D( ^2 2*Sl) , (15) 
k ^ k ^ 

with the shift operator 0(E^4) = ex P fefc(W w iO(' ! >*: ( 4 _ Skdkj), bath hamiltonian iJ ba th = J2k ^ka\a k , 
and the effective interaction 

H BS = -T,&S k Sl (16) 

k ^ 

It is very instructive to calculate this interaction explicitly for the example at hand. To that end we revert momentarily 
to the standard expression of the electric field, ea. (|21f) . and allow for arbitrary orientation of the DQDs. If we denote 
the orientations by unit vectors Hi and observe that the index k in ea. (|14f> stands for wave vector and polarization 

(k, a), we have Ska = S<=i i 5ka <J 2* elkri with coupling constants 



(i) , / go|kj 

■ 9k « = di ' ekQ V^' (17) 

with dj = (ed/2)uj. The evaluation of the sum over all modes in ea. l|16|l leads, in the limit of continuous k and infinite 
cut-off frequency, to 

d 1 -d 2 -3(d 1 -f)(d 2 -f) 
H cS = ^-^3 a zl a z2 , (18) 

i.e. a dipole-dipole interaction. This explains the 1/r 3 dependence of the phase <p~(t, to) and its proportionality to 
t for large times. Of course, the electric field can not mediate an instantaneous interaction, as ea. l|18fl might suggest. 
The full time dependence must take into account also the free bath hamiltonian, and iibath + H c ^ together lead indeed 
to the correct retardation behavior. This is seen from ea. l|12|) when we take the limit w max — * oo. The first oscillating 
term sin(y max to) arises from the sharp cut-off and will average to zero for a smoother cut-off. The last term vanishes 
for y m ax oo, and the remaining sin-integral functions conspire to a Heaviside theta-function on the light-cone, 
n6(t/to — 1), with the consequence that no entanglement can be created faster than the speed of light with this effect. 
The matrix C in eq.@ is seen to arise from the difference of eigenvalues of a z ia z2 defined by the indices of the density 
matrix. 

One might be tempted to think at this point that the entanglement generation is trivial in the sense that the 
hamiltonian for the two DQDs and the heat bath just amounts to a fancy reformulation of the dipole-dipole interaction 
between the DQDs initially assumed non-interacting. This is, however, not the case: First of all, as can be seen from 
eq. (|15f) the effective interaction is not the only term that determines the dynamics of the entanglement generation. 
Rather it is supplemented by two shifts in the harmonic oscillators, which depend on the state of the two DQDs. The 
reduced overlaps of the harmonic oscillators are responsible for the decoherence quantified by the functions fi(t,to) 
and /2(t, to). Thus, there needs to be a balance between the effective interaction induced and the decoherence due to 
state dependent coupling to the heat bath, and this balance is made quantitative by the functions fi(t,t ), ./^(Mo), 
and <P-(t, to). The effect of decoherence (and possibly the retardation) would be overlooked, if one started directly 
with a hamiltonian containing the dipole-dipole interaction. 

Secondly, the heat bath is initially in thermal equilibrium, and the thermal noise reduces the amount of entanglement 
additionally compared to the T — case. So the nature of the environment as a heat bath is important, and not just 
the fact that it induces a well-known interaction. 

One might wonder whether one can speed up the entanglement generation by reservoir engineering. One obvious 
attempt would be to selectively supress the coupling to one of the two heat baths, say the sin-waves. In fact, the 
calculation shows that for a complete supression of the coupling to the sin-waves the entanglement becomes basically 
independent of R and entanglement would be created quasi instantaneously, FigEJ The difference arises from the fact 
that ip- = cp i — <f2, a phase difference accumulated from the couplings to the cos- and sin- waves, is replaced by a 
single phase from the cos-waves alone. The latter contains a large distance independent term which normally cancels 
a corresponding term from the sin-waves. Without the sin-waves (ip2 — 0) this term leads to a very rapid, distance 
independent growth of (p- and therefore to the quasi instantaneous entanglement creation. We conclude that the two 
independent heat baths counteract each other when it comes to entanglement creation, and this leads to much slower 
entanglement creation based on the small remaining, distance dependent phase accumulation. 
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FIG. 3: Entanglement of formation E for a hypothetical coupling to the cos-waves only. Entanglement production is basically 
independent of R in the interval shown, and thus possible quasi-instantaneously (t < to). Same parameters and grey-scale code 
as in FigQ 



While a comlctc suppression of the sin-waves (or the coupling to them) might be illusive, even a very small suppres- 
sion (as one might imagine in a cavity by using a thin, uncharged wire) would lead quickly to entanglement creation 
faster than the speed of light. Fig0] shows the entanglement for to = 10 6 r as function of 7 and log 10 (t/r), where 
< 7 < 1 measures the relative coupling strength to the sin-waves (/^(ijto) ~~ * 7/2(^,^0) and <f2(t,to) — > 7(^2(^,^0))- 
The smaller 7, the faster the entanglement arises, but up to 7 very close to unity, almost perfect entanglement is 
created for t < to- 

However, causality leads to strong restrictions of what should be possible in this respect: also in a classical field 
theory the retardation of the dipole-dipole interaction would be modified, and a electromagnetic interaction spreading 
faster than the speed of light is certainly not possible, regardless of whatever arrangement of conductors one might 
come up with in order to supress certain modes. 

In (la , ITiil ] it was proposed that two atoms can get entangled by coupling them for a time t < to to the vacuum 
of a massless scalar relativistic field (thus explicitly avoiding the effects of any effective induced interaction), and 
an experiment was proposed in 20] to demonstrate the corresponding effect in ion-traps, namely the entangling of 
the internal degrees of freedom of two ions in an ion-string faster than the speed of sound in the ion-string. The 
vacuum case can be retrieved from the calculation presented here by setting T — 0, and replacing the coth function 
in ea. H10|) by unity. Since r diverges then, it is more reasonable to use directly to as time scale. This is achieved 
formally by substituting y — zr/to in eas. Hl(Jlll|) whereupon the functions depend only on t/to up to an additional 
prefactor r 2 jt\. The latter combines with the prefactor A and reads then aod 2 / (-/rcgi 2 ,)- The cut-off y m ax is replaced 
by z max = Lo max to- However, nothing changes about the fact that the function is proportional to a Heaviside-theta 
function centered on the light-cone. Thus, no phase accumulation is possible for t < t , unless one uses a finite cut-off. 
In this case the Heaviside theta function gets "softened", which is physically plausible as now the heat-bath does 
not contain sufficiently small wave-lengths anymore to precisely locate Bob's DQD. But even if in such a case one 
was to operate slightly before the light-cone, t < to, one would still have to overcome the smallness of the prefactor 
aod 2 1 (ncoto) ■ This seems to be excluded for atoms, as the distance between Alice and Bob would have to become 
smaller than the length attached to the atomic dipole moment d, in which case the whole dipole approximation breaks 
down. In DQDs one might hope to create states |0) and |1) with large amounts of charge, thus increasing artificially 
the dipole moments, but even for ten excess electrons in one well the DQDs would have to be so close that the dipole 
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FIG. 4: Entanglement of formation for to = 10 6 r as function of 7 and log 10 (i/r). The parameter 7 measures the coupling 
strength to the sin-waves (0 no coupling, 1 full coupling). Already values of 7 only slightly smaller than 1 would lead to 
entanglement creation before the light cone (t < to). Same parameters and grey-scale code as in FiglTI 



approximation becomes doubtful. Thus, within the above theoretical framework of fixed DQDs with two degenerate 
energy levels dipole-coupled to BBR, there seems to be not much room for substantially entangling the atoms faster 
than the speed of light by coupling them for a time t < to to the BBR, even at T = 0. 



VI. CUT-OFF DEPENDENCE 



As the function 93 _ remains finite for uj max — > 00, a different cut-off function has little influence on ip_ (besides the 
eventual removal of the oscillating terms mentioned earlier) and therefore on the speed of entanglement generation. 
However, since /i^ diverge for y max — > 00, the cut-off function can change the maximum amount of entanglement. 
Let us suppose that the cut-off function to be used in eq.JTUJ is equal to unity for y < y max and equal to a function 
C(y) for y > y ma x- For y max 3> 1, the coth(y/2) can be replaced by unity without impacting the following scaling 
arguments. Also, the strongest diverging i _ dependent part is given by 

dyC(y) — , (19) 

IX y 

which even for C(y) = 1 remains finite. The only remaining question is then, if a different cut-off function changes 
significantly the i-dependent behavior of /i,2, and in particular, if for t — > 00 /12 always remain finite. To answer 
this question observe that for < C(y) < 1, 

< - / dy C{y)y{\ - cos(yt)) < - / dy C(y)y . (20) 



"Umax - 

Thus, if the cut-off function decays fast enough to make the integral finite, it will be finite for all times. More 
specifically, if C(y) oc l/y p , p needs to be larger than two, and the upper bound in ij2"UI) is then of the order y 2 ^^, 
which is completely negligible compared to the dominating y^ax behavior from y < y max . Therefore, the results 
obtained, in particular the entanglement for large t, are robust against a change of the cut off function. 
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VII. CONCLUSIONS 



The CMB or the black body radiation (BBR) in a box-shaped cavity can be regarded as two independent heat 
baths, which couple differently to two spatially separated DQDs. Two spatially separated, non-interacting quantum 
systems can get entangled by interacting with the CMB or the same BBR. However, the effects of both heat baths 
cancel to large extent. As a result, entanglement is created only very slowly, with a first maximum entanglement 
arriving at a time ti that scales like the 3rd power of the distance R between the quantum systems. The effect can be 
understood as originating from an effectively induced dipole-dipole interaction, and consequently is retarded by the 
propagation of a light signal. 

Acknowledgment: I thank Jaewan Kim and Sangchul Oh for useful discussions, and Benni Reznik for useful corre- 
spondance. 



VIII. APPENDIX: THE ELECTRO-MAGNETIC FIELD 



I derive here the expression for the electric field operator, eq.(J2Jl. Starting point is the standard representation of 
the quantum operator of the electric field for periodic boundary conditions |34| , 



E ( r ' f ) = ^2 [a^ a e 1 



-i(w fc i-kr) , t „i(u)fct-kr) 
~r u ka e 



(21) 



with £k — y/hujk/(2eoV), unit polarization vectors e-^ a (a — 1,2), and the creation and annihilation operators a kcc 
and Ok Q j respectively. We express these operators in terms of canonical coordinate and momentum operators q^ a and 
Pkcn respectively, to obtain 



E(r, t) = V* e ka J_ (oj k qka cos(kr - Lu k t) - p k a sin(kr - Lu k t)) 



(22) 



We then split the set of modes into two sets, one with k x > 0, the other with k x < 0, and introduce Q^ a = 
(<?ka i <Z-ka)/v / 2, P k = (pket ±P-k«)/v / 2- Changing the summation variable from k — * — k in the k x < part, and 
with eki = — e_ki, ek2 = £-k2, si = — , S2 = +, Eq.JSlJ can be rewritten 

(w fc Q k ° coskr - P^ a sinkr) coscj fe t + (uJkQ^~ a sinkr + coskr) smto k t 



k,a 



where the prime at the sum denotes summation over modes with k x > only. Note that hereby the total number of 
modes is kept unchanged. We finally perform a time dependent canonical transformation, 



Qkii(i) = Qki coswfct + — P kl sinw fc i 

<3ki2(£) = QkiSinwfet — Pjtcosw/fet, 

u k 



(23) 
(24) 



and the same set of equations, but with Q kl — ► Qk2' ^ki ~~ * ^k2 ^ or ^ ne second polarization direction, Qk2i, Qk22- 
This leads to cq.@. The advantage of this representation is that all modes are coupled via their canonical position 
to the DQDs. 
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